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Two-fluid matter-quintessence FLRW models: energy transfer 
and the equation of state of the universe 
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Abstract. Recent observations support the view that the universe is described by a FLRW model with ~ 0.3, 
~ 0.7, and w < —1/3 at the present epoch. There are several theoretical suggestions for the cosmological A 
component and for the particular form of the energy transfer between this dark energy and matter. This gives a 
strong motive for a systematic study of general properties of two-fluid FLRW models. We consider a combination 
of one perfect fluid, which is quintessence with negative pressure (pq = weQ ), and another perfect fluid, which is 
a mixture of radiation and/or matter components with positive pressure (p = /3e m ), which define the associated 
one-fluid model (p — je). We introduce a useful classification which contains 4 classes of models defined by the 
presence or absence of energy transfer and by the stationarity (w = const, and (3 = const.) or/and non stationarity 
(w or (3 time dependent) of the equations of state. It is shown that, for given w and j3, the energy transfer defines 
7 and, therefore, the total gravitating mass and dynamics of the model. We study important examples of two-fluid 
FLRW models within the new classification. The behaviour of the energy content, gravitating mass, pressure, and 
the energy transfer are given as functions of the scale factor. We point out three characteristic scales, oe, a-p and 
a,M, which separate periods of time in which quintessence energy, pressure and gravitating mass dominate. Each 
sequence of the scales defines one of 6 evolution types. 
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1. Introduction 

A number of recent observations reveal the cosmolog- 
ical A component. Type la Supernovae (Riess et al. 

|200l| ) and the 



199?, Perlmutter et al. 1999, Riess et al 



Boomerang, Maxima and Dasi measurements of the total 
density parameter Q via the first acoustic peak location 
in the angular power spectrum of the CBR (de Bernardis 
et al. |2000| , Balbi et al. |2000| , Jaffe et al. |2000[ ) show that 
n = fij£ + Q° A = 1 ± 0.05 with Q° A « 0.7. Existing cosmo- 
logical data allow a wide range for the equation of state 
coefficient w fro m —1 to —1/3 (Pe rlmut ter et al 
Podariu & Ratra |2000j Wang et al. |2000|) . 



1999a 



Also, the smooth Hubble flow around our Local Group, 
inside a highly lumpy matter distribution, suggests still 
other evidence for a dominating A component and its 
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var iation with ti me (C hernin et al. 2000]; Baryshev et 
al. |200l|; Chernin |200l[ Klypin et al. |2001t Axenides & 



Perivolaropoulos |2002| ) 



There are several theoretical models for A-like cos- 
mological components of the universe with positive en- 
ergy density and negative pressure, including vacuum 
with a constant A, decaying A, and variable equat ion of 
state w(t) (Pee bles & Ratra |98§ Lima & Maia |l994 ; 



Wcttcrich 1995; Ferreira fc Joyce 1997; Caldwell et al 



1998; Stcinhardt et al. 



1998; Zlatev et al. 1999; Bahcall et 



al. 1996 ; Mbonye 2002 ). We use the term "quintessence" 
for any kind of substance having the equation of state 
Pq = w eq with — 1 < w < 0, which may be time variable. 

Thus observations and theory give strong motiva- 
tion to study general properties of two-fluid Friedmann- 
Lemaitre-Robertson-Walker (FRLW) models in which 
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the A component dominates at late epochs and there is 
energy transfer between A and matter components. 

The energy transfer between dust-like matter and radi- 
ation was first studied by Davidson (1962). Some examples 
of the physics producing the energy transfer were given in 
Sistero Ql97ip and Mcintosh Q1967] , |196^ ). The first exact 
solution of the equation of motion for a dust+radiation 
model with no energy transfer was obtained by Chernin 
(1965), who applied it to the case where the radiation com- 
ponent is the cosmic background radiation or neutrinos. 

In the present paper we give a systematic presentation 
of the properties of two- fluid cosmological models, with 
and without energy transfer, in the frame of a new clas- 
sification which naturally arises in the two-fluid problem. 
In section 2 we briefly summarize the standard two-fluid 
FLRW model. In section 3 we consider general properties 
of two-fluid models with matter and quintessence. Section 
4 includes three examples of models from different classes 
of our classification. Section 5 contains the conclusions. 



2. The two-fluid model: a summary 

The derivation of the FLRW equations contains the follow- 
ing basic elements: 1) The Einstein equations (notations 
from Landau & Lifshitz 1971) 

w k -\glK = ^n (i) 



s s 2 



kc 2 



(2) 



2) The Bianchi identity in the form: 

3) The RW line element in spherical comoving space co- 
ordinates (x, 9, 4>) and synchronous time t: 

ds 2 = c 2 dt 2 - S 2 (t)d x 2 - S 2 (t)I 2 ( x )(d6 2 + sin 2 ddcj) 2 ) (3) 

4) The total energy momentum tensor in comoving coor- 
dinates (ordinary matter, vacuum, quintessence): 



T, 



k - diag(e, -p, -p, -p). 



(4) 

Here Ik(x) = ( s * n (x)i Xi sinh(\)) for the curvature con- 
stant k = (+1,0, —1), respectively. S(t) is the scale factor, 
e = pc 2 is the energy density, and p is the pressure. 

The proper metric distance r of a particle with the 
dimensionless comoving coordinate x from the observer is 

r(t,x) = S(t)- X . (5) 

Here physical dimensions [r] = [S] = cm. The volume 
element in metric (||) dV — S 3 1 2 (\) sin 2 (9) d\ d9 dcj). 



2.1. The FLRW one-fluid model 

For the one-fluid model the (0,0) and (1,1) components of 
Eq.(|l|) give the FLRW equations 

c-2 

(6) 



S 2 ' 'S 2 



8ttG 

•— p. 



S S 2 S 2 <? 
where ' = d/dt. Eq.(||) implies: 



(7) 



(8) 



S e+p 

The equation of state for the one-fluid model is 

p = je. (9) 

The characteristic energy content e of the sphere with the 
radius equal to the proper metric distance r = S(t)x is: 

(10) 



e(r) = / edV ^ineS 3 a k ( X ), 



where for different values of k = +1, 0, —1, <Jk(x) = \ — 
4, _ 2^ respectively. 

The gravitating mass M g is defined by the equation of 
motion, which follows from Eqs.(||) - (0): 



-G 



M g {r) 



where 



M 9(r) = £ -^ E I = 4(37 + 1). 



(11) 



(12) 



This definition corresponds to the general expression 
for the gravitating ma ss of a self-gravitating body as given 
by Landau & Lifshitz |l97l|, §100). 



2.2. Time dependence of energy 

There are two factors which cause a time dependence of 
energy of each particular fluid and, hence, of an associ- 
ated (as termed by Coley & Wainwright 1992| ) one-fluid 
model. The first cause is the expansion of the universe. 
The second is when one fluid is converted into another. 

1931 ; Davison 1962c Harrison 



It is know n (Lc maitre 



1995| ; Peebles |l993j p. 139) that Eqs.(|) and (|7j) imply the 



non-conservation of the energy s of one-fluid FLRW mod- 
els filled by a non dust-like matter: 



de d , o, 
M +P dt {R) =°- 



(13) 



For every particular fluid the time dependence of the en- 
ergy is defined by its pressure. Eq.(|l3|) is a consequence 
of the Bianchi identity (||). Rewritten for each particular 
fluid, it has the form of the classical energy conservation 
law without the convertion of one-fluid to another: 



3c 2 



de Q 
dt 



0. 



ds r 



dt 



0. 



(14) 
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When the one-fluid converts into another (for instance, 
the matter into radiation in stars), Eqs.(p^) become 



1 

R 3 

1 

R 3 



dt 
dt 



»4< R3 > 



UQ, 



■ PmJt {R*) 



= 0, 



u m and uq are the rates of energy transfer per unit volume 
from the quintessence to matter and vice versa {[u m \ = 
erg/sec-cm 3 ). The terms pq-^(R 3 ) and p m -^(R 3 ) describe 
the time dependence of the energy due to space expansion. 

2.3. The two-fluid model for matter and quintessence 



In his pioneering work, Lemaitre ( |1933[ ) studied a model 
containing dust and the cosmological constant. The first 
application of the ener gy transfer for two-fluid cosmology 
was by D avidso n ( 1962 ). We also m ake us e of the results by 
Chernin ( |1965| ), Thomas & Schulz ( |2000D , Coley & Tupper 
( |l986a|) , Amendola & Tocchini-Valentini (|200l| ). 

We study two groups of substances: 1) a number of 
perfect fluids with positive pressure; 2) a number of flu- 
ids with positive energy and negative pressure. Each group 
plays the role of a one-fluid in two-fluid models. The equa- 
tions of state for the first group are: 



Pn 



(16) 



where p m = p m i+Pm2+— , £ m = s m i+e m2 + ..., < fa < 1 
and /3 0. The second group has: 

(17) 

, -1 < 



PQ = we Q , 

where p Q = p Q i 
uij < and w = 



+ PQ2 + 

^ < 0. 

PQ 

We define an important function 



£ Q — e Qi + e Q2 



a(a) 



£q{o) 



(18) 



which characterizes the relative contributions to the en- 
ergy density from quintessence and matter (see e.g. 
Wetterich 1995). The relative energy contribu tion w as 
studied in another way by Coley & Wainwright ( 1992 ). 

The two-fluid model is defined on two levels. Firstly, 
on a level of partial one-fluid models, with two kinds of 
interactions, 1) through the common gravitational field 
and 2) by the energy transfer from one-fluid to the other. 
Secondly, on the level of the associated one-fluid model. 

For a complete description of a two-fluid model one 
has to study both levels simultaneously. We assume state 
equations of unique form for partial models and for the as- 
sociated one-fluid model. The essential difference is that 
the equation of state of the associated model is not pro- 
duced by physical particles and their interactions. 



2.3.1. The associated one-fluid model 

We specify the partial energies, pressures and the corre- 
sponding equations of state: 



£Q = PQC ■ 
PQ = we Q , 
(15) P = 7£- 



Pm pEffi , 



- £q + e m , 

PQ+Pm, 



(19) 



Thus the equation of state for the associated one-fluid is: 
' (20) 



p wa 
^ £ a 



1 



For the energy ratio a > one finds that w < 7 < (3. 
The associated one-fluid description of the two-fluid model 
includes Eqs.(||) - (0) and Eq.(pii|). The three functions of 
the scale factor a w(a), /3(a) and a(a) are assumed given. 

2.3.2. The two-fluid model 

The energy-momentum tensor of the associated one-fluid 
model is the sum of the partial tensors of matter and 
quintessence. According to Eq.(0) we have 



rpl rji I 

k-.i (m)k;i 

Eq.(||) implies: 

£tot 



T (Q)k-,i - °; 



(21) 



3§ = - 

<-> £tot + Ptot 



£tot = SQ + £m, Ptot = PQ +Pm- (22) 

There are two subcases within the two-fluid model de- 
scription. In the first the energy-momentum tensors of the 
partial fluids are separately conserved: 



T, 



= 



and 



T, 



(Q)k-i — 0' 



(23) 



(m)k\i 

so that uq = and u m = 0. In the second case we allow 
the presence of the energy transfer, when uq ^ and 
u rn ^ but uq + u m = which is equivalent to Eqs.(p2|). 

Each fluid has an equation of state in the form of 
Eq.([)|) with different coefficients: /3(a) for the matter com- 
ponent and w(a) for quintessence (see Eqs.(|l9|) ). 

For two-fluid FLRW models the characteristic energy 
content within the sphere of a fixed comoving radius x is 



4ireS 3 a k { X ) = 4tt (a + 1) e m S 3 a k ( X ) 



(24) 



The gravitating mass M g is the sum of the partial gravi- 
tating masses of two-fluids: 



M g (r) = M™ + 



M Q 

g 



[(1 + 3w)e Q + (1 + 3/3K 



dV 



(3 7 +l). (25) 
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2.4. Dimensionless equations 

Now we restate the models in terms of dimensionless quan- 
tities. We introduce the following characteristic values: 



1 _ 3ff 2 

-fig 8 7T (_r 



(26) 



4 7T 3 C 2 

M °-T P0 ^ = 2H CJ G' 
and dimensionless variables: 



e = Poc , 



Eqs. (§) - ® and @ become 



dak 
2- + + = -3P, 



P = 7 £, 



(27) 

(28) 

(29) 
(30) 



Table 1. Classification of two-fluid FLRW models. 
ET means "energy transfer"; NET is "no energy trans- 
fer" ; SES is "stationary equation of state" ; NSES is "non 
stationary equation of state" . 





u Q = u m = o 


U Q = ~U m + 


both w, (3 
axe constant 


NET-SES 


ET-SES 


at least one of w, @ 
is not constant 


NET-NSES 


ET-NSES 



a(r) 


S(t) 
lo ' 




P_ 

Po 1 








p(a) 


M g (t) 
M 


= 3<7fc(x)(37 - 


- l)a 3 £, 








£ = 


s 

£o' 


c e Q 


£"m — ) 

£o 


£-- 


= £q* 




P = 


P_ 

£o' 


P _PQ 


p _ P" 1 

r m. — , 


P 




Pm 


E = 


e 

eo' 


eo 


„ e m 
eo 


E 


= Eq 


E E n 


U = 


u 

— to, 

e 


Uq = — t Q , 

£o 


U m = — to, 


U 


= u Q 


f U m 



3. General properties of two-fluid FLRW models 
with matter and quintessence 

3.1. The classification 

A two-fluid model with the vacuum or quintessence to- 
gether with some mixture of matter is very different from 
the two-fluid model with "ordinary" matter: negative pres- 
sure and gravitating mass components give rise to a new 
behaviour of the total pressure and total gravitating mass. 
In order to facilitate their study, we divide all two-fluid 
models into four classes according to two independent 
properties: 1) both fluids have a stationary (SES) or at 
least one has a non stationary (NSES) equation of state 
and 2) the presence (ET, energy transfer) or absence 
(NET, no energy transfer) of an energy transfer between 
the components (see also Table 1). 

Furthermore, we separate three different kinds of two- 
fluid models depending on the behaviour of the function 
a(a) = £Q/£ m : the coherent model, when 



a(a) — const, 

the asymptotically coherent model, when 
lim a' (a) = 



(34) 



(35) 



and the non coherent model, when neither Eq.(34) nor 
Eq.(f35|) are valid. 



Uq 



Un 



d 

' d_ 



£ Q a 3 ) 



( &Ti 



;) d 3 



Prr 



dT 



3.2. The general solution of two-fluid FLRW models: 
ET-NSES models 



Eq.(|lT|) for the scale factor now has the form: 
a 2 d = — Oq n(a), 



(32) 



(33) 



where fio is the cosmological density parameter at the 
moment of the initial conditions for Eq. ( 33 ) . 



The ET-NSES class is the most general class of two-fluid 
models in our classification. We assume that all fluids (the 
associated and the two particular ones) have nonstation- 
ary equations of state. 

In this subsection we completely describe the two-fluid 
problem in terms of the associated one-fluid model and 
two particular fluid models simultaneously. This leads us 
to rewrite the FLRW equations in terms of the coefficients 
of the three equations of state: 7, (3 and w. 
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3.2.1. The input equations and their solution 

We describe the associated one- fluid model by Eqs.(j2£ 

e 



Another form of Eq. (|40|) is obtained by substituting 2 
from Eq. (|27j) . This shows how the energy transfer depends 
on the curvature (k ^ 0): 



3* 
a 



£ + P' 



P = 7 £, 7 = 7(a). 



(36) 



c a V a z 



a 
a 



(43) 



The solution of Eqs.pq), with the initial conditions stated 
for the present epoch, is: 



£ = £ (1) exp 



7+1 



dx 



(37) 



Using the definition of a, Eq.([18|), and its time 
derivate, a (7 — w) 2 = (7 — w) (3 + ((3 — 7) u> — ((3 — w) 7, 
Eq.(^o|) leads to the expression for the energy transfer: 



U m ("f,a) 



In this model the characteristic energy content E, pressure 
P and gravitating mass \x are: 



±3 



(/3-7)(7- w) 
(3 — w 



E= £a 6 



The scale factor satisfies Eq.(|33|). 7(a) is given by 
Eq.(po|) in case of known w(a) and /3(a). 



£7, Li = 3a k ( X )£ a 3 (3 7 + l). (38) (0 - w) 



(7 - w )p + {p- 7 ) v/-(p- w) i 



•(44) 



The two particular fluids are described by Eqs.(31) 
( |32| ) and the equations of state: 



P Q = w(a) £q, P rn = 13(a) £„ 



(39) 



Calculating a 3 (u m — ^rj and keeping in mind that 
Uq + U m = 0, we find the equation 



U m a + 1 3 
^ a 

-C-m O 



3-(3-w) . 

a a 



(40) 



Eq. ( f40| ) leads to a general property of the model: the model 
can be coherent (a — constant) only if there is energy 
transfer. Eq.(f4(i|) defines the critical energy transfer U 
which keeps a constant (the coherent solution): 



E 
m ' 



U*=3~£(J3- 
a 



uu) 



a 



(a + 1) 2 



u ™(w,/3,7, a )- 



(41) 



separates the behaviour of the two-fluid model with 
different signs of a(a): it follows from Eq.(|40|) that a(a) > 
requires U m < U^;; a(a) < requires U m > U^. In 
other words, the rate increases, if U m < and 
decreases if U m > U^. 

In a number of works the energy transfer has been cal- 
culated for a given model of interaction between dust and 
radiation, i.e. for given w(a) and (3(a) (see, for instance, 
Mcintosh |1967| and |g6§ Sistero |l97l| ; Coley & Tupper 
). Now we know that energy transfer is 



1985 and 



1986 



necessarily required for the coherence. 

At the same time Eq. (pp|) gives the condition for no 
energy transfer, Uq — U m = 0, in the form: 



a(l) exp 3 



-dx 



(42) 



Any other function a(a) leads to non zero energy transfer. 



£ (7, «) \ £{l,a) o' 



where ' = ^ and £(7,0) is defined by Eq. (|37|) . The sign 
'+' corresponds to expansion while the sign ' corre- 
sponds to collapse. So, we found that the energy transfer 
U m depends on 7, f3, w, their derivates and, also, on the 
scale factor a (as on an independent variable) and cur- 
vature. For a more realistic situation with given k, (3(a) 
and w(a), Eq.(f44|) demonstrates that 7(a) is completely 
defined by U m (a). 

Eqs.(fiol), (p^ ) and ( ^4| ) represent in different forms how 
a convertion of one-fluid into another (i.e. U m ^ 0) influ- 
ences the properties of particular fields (its equations of 
state, and its energy transfer rates) and the expansion of 
space (i.e. a(t)). Eq.([44|) states the dependence of the en- 
ergy transfer on 7 and scale factor as an independent vari- 
able. Excluding 7 from Eqs.flU), @ and @ we find that 
energy transfer completely defines the gravitatng mass as 
a function of the scale factor and, so, the dynamics of 
the universe (see Eq.([j3|)). Therefore, the total gravitating 
mass and, consequently, the dynamics of the expansion, is 
determined by the energy transfer. 

For the popular flat model Eq.(Eol) becomes: 



tt ( a + 1 ) 2 e m i a \ 



a 



(45) 



The general (ET-NSES) solution for particular fluids 
is represented through the function a(a) and two effective 
coefficients of the equation of state (3(a) and w(a): 



Eq 
Po 



_£a_ 

oTT' 

£ Q a\ 
w £q, 



f 



£ 



a + 1' 

fJ-Q = 3 (Tfe(x) (3/3 + 1) a 3 £q, 
fi m = 3 (Tfe(x) (3 w + 1) a 3 £ w . 



(46) 



6 



Please give a shorter version with: \authorrunning and/or \titilerunning prior to \maketitle 



a(a, 



ed ) 



(47) 



Depending on the initial conditions and the class of a 
model, there is one (or, possibly, more) moments of time 
when the one-fluid model is effectively dust, i.e. j(a e d) = 
0; a e d is the "effective dust scale" . The general expression 
for this scale is: 

A 

w 

We now show the complicated character of the evolu- 
tion of the two- fluid model from initial conditions, defined 
at a m i n , to an asymptotical state with a — > oo. To do this 
we study the relative contributions of energy, pressure and 
gravitating mass of each component to the total energy, 
total pressure and total gravitating mass. Eqs.@) give 
the criteria for the relative contributions in the form: 



Table 2. Detailed description of the EPM evolution type 



a 

V = 

M = 
where 
ap = 

a M = 



£m 

\Pq\ 

Pm 

ImqI 



a 

a-p ' 
a 



(48) 



otM{a) 



Ha)|' 
3 /3(a) 



1 



., , , , ,,- (49) 
\3w{a) + 1 

Three characteristic scales as, a-p and aM appear so that 
a = as implies a = 1 , 
a = a-p implies V — 1, ap = a ec i, 
a = aM implies M. = 1. (50) 
These are the solutions of the equations: 
a(a E ) = 1, 

\w(a P )\ 



(3{a P ) 
\3w(a M ) + 1| 



■ a op 



1, 



"(O-A/f) = 1- 



(51) 



3P(a M ) + l 

/3(a), w(a) and a(a) are model-dependent, so there is no 
possibility to study the general case. 

Eq.([h)|) shows that the energy transfer is required for 
the coherence (a — const). Now we give two other ex- 
amples of how to use E q.([40| ) . Substituting a = ap V 
and a — aM M into Eq.(40) we find for V = const and 
M. = const two critical values for the energy transfer: 

Va-p 



U p = 



3 - (/3 - w) 

a ar 



{Vap + 1) 2 ' 



and 
U M = 



3^(/3 
a 



Un 



w) 

= U p gives P Q 



d In aM 



(It 



£ 



Mum 



(M a M + 1) 2 



(52) 



(53) 



a-min < a < as 


£-m > £q 
Pm > \Pq\ 

MM > AJQ 


ae < a < a-p 


Em < £q 
Pm > \Pq\ 

UM > \fMQ | 


ap < a < a M 


£m < £q 
Pm < \Pq\ 

hm > ImqI 


aM. < a < oo 


£m < £q 

Pm < \Pq\ 

UM < 



3.2.2. The evolution types 

In this paper we study only the models which are fully 
dominated by the matter-component in the very early 
epoch (i.e. a < 1, V < 1 and £ < 1 ), while fully domi- 
nated by quintessence at the limit of large time (i.e. a > 1, 
V > 1 and £ > 1 ). Under this assumption the character- 
istic scales qe, ap and a M allow 6 types of nonequalities: 



a £ < ap < a M , 

as < aM < ap, 

ap < a M < a £ , 

ap < a £ < a M , 

aM < a £ < ap, 

a M < ap < a £ , 



EPM 
EMP 
PME 
PEM 
MEP 
MPE. 



(54) 



P m and U„, = U 



gives \i Q ~ ^„ 



Each system of nonequalities corresponds to an "evolu- 
tion type" , defined by a unique sequence of the scales a#, 
ap, and a M , which represents the sequence of time peri- 
ods when the quintessence-component dominates in total 
energy, pressure and gravitating mass, respectively. The 
ratio between two energies Eq (actually E\) and E m was 
first considered by Chernin et al. ( 2000| ) to make a divi- 
sion line between the growth and suppression of structure 
formation. In Table || the EPM evolution type is described 
in detail. Eq.(pp|) leads to a dependence of a on U m for 
given w(a) and /3(a). Eqs.(^|) - ( |5l| ) lead to a dependence 
of the evolution type on a for given w(a) and /3(a). So, we 
resume that for given w(a) and /3(a) the evolution type 
is connected with the energy transfer. On the other hand, 
energy transfer defines the gravitating mass. So, there ex- 
ists a correspondence between the dynamics of a model, its 
evolution type and energy transfer. 

Below we study how quintessence influences the behav- 
ior of the total gravitating mass, pressure, and the effective 
equation of state of the associated one-fluid model and we 
calculate the evolution types for a few examples. 
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4. Classifying two-fluid FLRW models: examples 

4.1. An example of NETS ES models 

This class of two-fluid models has been much studied be- 
cause the input equations are simple. In this subsection 
we study the evolution types for w — — 1, —2/3, —1/3 and 
for (3 = 0, 1/3, 2/3, 1. 

4.1.1. Input equations and their solutions 

The input equations are given by Eq.(E3) for every fluid: 



£, 



Q 



£q + Pq 
<~ in 



3 _ = _ 

& £m "T" Pfn 

u Q = 0, U m = 0. 
Eq.(f40|) is written in the form: 

3-(l3-w)--=0 
a a 

and has a solution 



, Pq — w £q, w = const., 
■, P m = (3 £ m , (3 = const., 



(55) 
(56) 
(57) 



The solution of Eqs.(|55|) is well known. Energy density 
£(a), characteristic energy content E(a), pressure P(a) 
and gravitation mass u(a) have the form: 



E Q (a) 



£q0) 

a 3(w+l) ■ 



7 3 w 



£ m (a) - 
E m (a) 



£m(l) 

a 3(/3+l) ' 

£m(l) 

n.3/3 ' 



P Q (a) = w 



gq(l) 

a 3(w+l) ' 

fj, Q (a) = 3o*(x)£q(1) 
^rn(a) = 3a k (x)£m0-) 



^W^Jwi)' (58) 

3w + 1 



a 3w 

3(3+1 



,3/3 



. We comment on two models in this class. First, (w = 
— 1/3, (3 = 0) has a constant total gravitating mass and 
total scale-dependent energy content. Second, (w = —1, 
(3 = 0) produces a total constant pressure and non con- 
stant energy content and gravitating mass (see Table g). 

4.1.2. The condition of coherence 

Assuming coherence (see Eq.(p4h) Eq.(p7|) leads to 



([3-w)= 0, 



(59) 



which for the non-static solution reduces to the simplest 
case of one-fluid dust model: 



Assuming asymptotic coherence (Eq.(|35|)) Eq.(p7[) gives 

ft ' = 3a^ = 3a(l)(/3-«j)a 3 (' 3 -'")- 1 , (61) 
a 

so, lim^oo a' = only for (3 — w < |. 

4.1.3. Asymptotic behaviour and evolution type 

The effective coefficient 7(a) of the associated one- fluid de- 
pends on time. Substituting ( |57| ) into ( pfj| ) one finds 7(a): 



7(a) = (3 



(3 — w 
~T~ - 



1 



(62) 



which has the following asymptotical properties: 
lim 7(a) = (3, lim 7(a) = w, (63) 

a — >0 a — >oo 

so, any matter- quintessence NET-SES model is effectively 
a matter model at an earlier time and quintessence at the 
limit of large time. 

There is a unique moment of time when the total pres- 
sure is equal to zero and the associated one-fluid model is 
dust-like. The corresponding effective dust scale a e a is: 



a e d 



-w B 
~ A 



(64) 



To find the evolution types of NET-SES models, we 
use the function a to calculate two functions V and M: 



a(a) 

V(a) 

M{a) 



£m(l) 
I '"' 

13 w 



a(a), 
II 



3/3+1 



a(a). 



(65) 



The characteristic scales (see Eq . (pOQ) may be calculated 
analytically: 



Clp 



p £ m {iy 

WW). 

3(3+1 £JX) 
\3w + l\£ Q (l) 



(66) 



(3 = w = 0. 



(60) 



The numbers in Table 4 illustrate the complicated char- 
acter of the evolution from the matter-dominated model 
at the earlier time to the quintessence-dominated stage 
for large time. Different physical quantities pass through 
'boundaries' (a = 1; V = 1; M. = 1) between the matter- 
dominated and quintessence-dominated stages at different 
times. It is also seen that a ec i cannot completely charac- 
terize the evolution type and dynamics of the universe. 
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Table 3. Two examples of the NET-SES models: w = -1 
and = provide a constant pressure; w = —1/3 and 
= provide a constant pressure. 



w 





E 


P 


A* 


-1 





Aa 6 + B 


-A 




-1/3 





Aa + B 


-A/(3a 2 ) 


B 



Table 4. Evolution types for NET-SES models with u> = 
-1, -2/3, -1/3 and for = 0, 1/3, 2/3, 1. These numbers 
are obtained from Eqs.(|50|) and ( |5l| ) with £q(1) = 0.7 and 
£ m (l) = 0.3 We denote by (ME) the case of = as 
and so on. 





it) = —1 


it; = -2/3 


it; = -1/3 




ap = 0.754 


ap = 0.655 


ap = 0.429 




a P = 0.00 


a P = 0.00 


ap = 0.00 




a M = 0.598 


ajvt = 0.655 


aM = oo 




PME 


P(ME) 


PEM 




a ed — oo 


a ed — oo 


a ed = oo 


= 1/3 


ap = 0.809 


ap = 0.754 


ap = 0.655 




a P = 0.615 


ap = 0.598 


ap = 0.655 




aw = 0.809 


a M = 0.949 


aM = oo 




P(ME) 


PEM 


(PE)M 




o el j = 1.07 


a ed = 0.950 


a e d = 0.655 


= 2/3 


ap = 0.844 


ap = 0.809 


ap = 0.754 




a P = 0.778 


ap = 0.809 


ap = 0.949 




= 0.915 


um = 1-07 


ajvi = oo 




PEM 


(PE)M 


EPM 




o ed = 0.915 


a ed = 0.809 


a ed = 0.598 


0=1 


a E = 0.868 


ap = 0.844 


ap = 0.809 




a P = 0.868 


ap = 0.915 


a P = 1.07 




a.M = 0.975 


OjK = 1.11 


aM — oo 




(PE)M 


EPM 


EPM 




a ed = 0.868 


a ed = 0.778 


a ed = 0.615 



4.1.4. The time dependence of the gravitating mass 

The total gravitating mass in any NET-SES model has its 
extremum at 



30+l\^^> 
-w 3 w + 1 



(67) 



As a* should be real, so —1/3 < w < 0. The requirement 
of positive //'(a*) > is reduced to > u>, which is always 
valid. So, in NET-SES models the gravitating mass has a 
minimum for — 1/3 < w < and there are no extremums 
for -1 < w < -1/3. 

4.2. A NET-NSES example 

Here we generalize the known solution for mixed radiation 
and dust, first obtained by Chernin ( |1965 ). We use it as 
the matter component of a two-fluid. For quintessence we 



assume w — const. A new property, produced by the com- 
plicated matter-component is that the effective equation 
of state is non-stationary. So, this model is NET-NSES. 

4.2.1. The input equations and their solutions 

The input equations are given by Eqs.(p9|) for each fluid: 
to 



3^ 
a 



f — 



£q + Pq 

Sr 

Er + Pr 
_£d 

£d 
£r + £d 



P Q = w £ Q , 
P _ £ R 

Pd = 0, 

Pm =P D +Pr = P J 



w — const. 



U Q =0, Ur = 0, U R = 0. 



(68) 



The solutions of Eqs.(p8[) are well-known. The energy den- 
sity £{a), characteristic energy content E(a), pressure 
P(a) and gravitating mass fJ.(a) have for each component 
the form: 



£q = 

Eq = 



gq(l) 

a 3(w+l) ■ 



1 3 w 



Pq = w 



£qW 



a 3(w + l) ' 
3w + 1 

UQ - t Q {l) fl3m , 

For a and we find: 
a(a) = 

0(a) = 



, £ D (i) f 

CD — 3 , OR 



Ed — £n(l), Er 
Pd - 0, P R - 

md = £d(i), 



a 4 ' 

£rW 
i £r(i) 

3 a 4 ' 



S D (l)a + £ R (l) 
V m _ £r(1) 



E m £ D {l)a + £ R {iy 
Eq.© has the form: 

3 - (0(a) - w) - - = 0. 
a a 

4.2.2. The conditions of coherence 



{69) 

(70) 
(71) 

(72) 



From Eqs.(|3J), (|70|) and ([72]) we find 

, 3C* a 1 - 3 ^ 

a = ; 

a Aa + B 

lim a' = -3w^a- {3w+1) = 0, for w > -- 

a— »oo A 3 

Q C f 1 

= _3w 34' for w = _ 3 

1 



(73) 



for w < --. (74) 
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Table 5. These numbers are obtained from Eqs.(50) and 
(fF§) with £q(1) = 0.7 and £ m (l) = 0.3; we choose here 
£ D = 0.05, S R = 0.25 



w = —1 


w = -2/3 


w = -1/3 


a E = 0.803 
a P = 0.587 
a M = 0.788 

PME 
a ed = 0.587 


a B = 0.743 
a P = 0.563 
a M = 0.921 

PEM 
a e d = 0.563 


a E = 0.634 
a P = 0.598 
aM — oo 
PEM 
a ed = 0.598 



4.2.3. The asymptotic behaviour and evolution types 
of a dust-radiation-quintessence model 

The effective coefficient 7(a) of the associated one-fluid 
depends on time. Substituting ( f7(i| ) into (p0|), one finds: 



7(a) = 



£ Q (l)w 



gfi(l) ffi 3io-l 



E Q (l)+£ D (l)a 3 <°+£ R (l)a 3w - 1 



(75) 



Because /3 — w > 0, the associated one-fluid model has the 
following asymptotical properties: 



liiWa) = -, 



lim 7(a) = w, 



(76) 



So the one-fluid model, associated with the NET-NSES 
two-fluid under consideration is effectively raditive at 
early times and quintessence- like at the limit of large time. 
The model is effectively dust on the (unique) scale 



a e d 



-3w£ Q (l) 



(77) 



The evolution types are characterized by the three func- 
tions a, V and M.: 



V = 



3 w 
Ml) 



(£ D (l)a + £ R (l)) a(o), 



£D{t)a + 2b R (l) 



In Table 5 three characteristic scales represent two evo- 
lution types for this model: PME and PEM. 

We note that the two-fluid description coincides with 
the one-fluid model only asymptotically, at a — > 00. 

4.2.4. The time dependence of the gravitating mass 

For — 1/3 < w < the total gravitating mass has a unique 
minimum at a scale a* 



(79) 



4.3. An ET-SES example 

Unlike the NET-SES models, ET-SES models allow co- 
herence. We inspect here an ET-SES example, with the 
two-fluids described by constant (3 and w. 

An important class of dark energy (DE) models is the 



so-called coupled quintessence (see e.g. Wetterich 1995 
Amendola 2000), where the total (matter + DE) en- 
ergy momentum tensor is conserved, whereas in ordinary 
quintessence the matter and DE are separately conserved. 
For coupled quintessence, a stationary DE model was re- 
cently found (Amendola & Tocchini-Valentini [2001 ) which 
predicts a coherent behaviour of matter and DE densi- 
ties at late cos mic epochs. S uch a relation was studied by 
Baryshev et al. ( 2000a 2001 ) in order to solve the problem 
of the low velocity dispersion in the local Hubble flow. 

The coherence (a = constant) implies 7 = constant 
and the model is greatly simplified. All these properties 
are produced by the energy transfer. 

4.3.1. Input equations and their solution 

The associated one-fluid model is described by 

£ 



3- = - 
a £ + P 



P = 7 £, 7 = const. 



which have the solution: 
£{1) 



£ = 
P = 



a 3 (-r+i)' 
£(1)7 



E 



a 3 M a 3 ^ 

For the quintessence component we find: 
£(l) a „ £(1) a 



g(l) 

£(1)(3 7 + 1) 



£q = 

Po = 



a 3( 7 +l) a+ y 

£(l)w a 



Eq = 



a 3 (7+i) a + r 

n Q = 3cr fc (x) £(1) ■ 
and for the matter-component: 
£(1) 1 



a 3 T a + 1' 



a 3 w + 1 



p — 

P — 



a 3 (7+D a + 1' Em 

£{i)P 1 

a 3 (7+1) a + 1' 

Mm = 3cr fc (x) £{1) ■ 



+ 1 a 3 7 



£{l) 1 



a 3 7 a + 1' 



1 3/3+1 



a + 1 a 3 7 
From Eq.(p0[) we find the energy transfer 

a £(1) a 
Um = 3 ^TWa^I)a {P - W) - 



(80) 



(81) 



(82) 



(83) 



(84) 
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4.3.2. The condition of coherence 

The coherence implicitly exists in this model. 

4.3.3. The asymptotic behaviour and evolution types 

The asymptotic behaviour strongly depends on the sign 
and value of 7, obtained at the end of the section (4.3.1). 
Three functions, generally defining the evolution type, 

£n 

a = -r 2 - = constant, 



V = —- a = constant, 
P 



M = 



|3w + l| 
3/3 + 1 



a = constant. 



(85) 



are constant in this model (a fixed ratio between the en- 
ergy, pressure and gravitating masses of the quintessence 
and matter components). What component dominates 
during all the evolution, depends on the value of a. 

Let us study now the case of the Amendola & Tocchini- 
Valentini ( 2001 ) solution, which is reached at late times: 
w = —0.7, j3 = 0, a = 7/3. The equation of state of the 
associated one-fluid model is defined by 



7 



w a + (3 
a + 1 

which gives 

6(1) 



£ 



49 

Too 



6(1) 



(86) 



O 153/100 fl 3/2 

E = £(l)a 147 / 100 «£(l)a 3 / 2 



49 

V = - — 



6(1) 



100 a 153 / 100 



1£(1) 

2 a 3 / 2 



47 

a = -6(1) a 1 

h 100 v ; 



1 



-§£(i) « 3/2 - 



147/100 

The energy transfer for the flat (A; = 0) model is: 

Urn. — 3 ■ 



(87) 



! g 3/2 (l) i g s 

(« +1) 2 fl 9/2(7+l) ^ P W> 



This gives for the Amendola & Tocchini-Valentini solution 
441 £ 3 / 2 (l) 2 £ 3 /2(l) 



U, 



1000 a 459 /2°° 



,2.3 



(89) 



has initiated many studies of FLRW models with an essen- 
tial A component at late epochs. Usually one has viewed A 
and matter as independent substances so that the energy- 
momentum tensors of the partial fluids are separately con- 
served (Eq.(p3|)). However, there are a number of sugges- 
tions in the litterature on the particular forms of energy 
transfer between dark energy and matter. This has mo- 
tivated us to consider on a phenomcnological level the 
general case when one-fluid converts into another and the 
equation of state for both components is non-stationary 
(Section 3). The properties of the model we gave in terms 
of the coefficients of the equation of state for two partial 
models (/? for a component with positive pressure and w 
for a negative pressure one) and a coefficient 7 for the as- 
sociated one-fluid. The energy transfer, when a one-fluid 
converts into another, was represented via these coeffi- 
cients and the scale factor (Eqs.(^) - ([45])). 

We have analyzed four classes of models defined by 
the presence or absence of the energy transfer and by the 
stationarity (w — const., and = const.) or/and non- 
stationarity (w or (3 time dependent) of the equations of 
state (see Tables 1, 4, 5). It was shown in sect. 3 that 

* for given w and /?, the energy transfer defines 7 and, 
therefore, the total gravitating mass and the dynamics 
of the model. 

* also, the model can be coherent only if there is energy 
transfer. 

The classification was illustrated with interesting exam- 
ples of two-fluid FLRW models in sections 4.1 - 4.3. 

From the behaviour of the energy content, gravitat- 
ing mass, and pressure as functions of the scale factor we 
have defined three characteristic scales, ole, a-p and aM- 
These separate time intervals when quintessence energy, 
pressure and gravitating mass were dominating (Eqs.(|48|) 
- (|5l]) ) . Any sequence of the scales defines one of 6 evolu- 
tion types of the model (Eqs.jS^)). There is a correspon- 
dence between the dynamics of a model, its evolution type 
and energy transfer. 
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Thus this model is asymptotically NET-SES. 



5. Conclusions 

The cosmological view that the universe is described by a 
FLRW model with fl a m w 0.3, n° A w 0.7, and w < -1/3 
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